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PREFACE 

Part of the contribution of the Ontario Minis tiry of 
tile Environment to the International Field Year for the 
Great Lakes (IFYGL) program was the estimation of ground- 
water inflow to Lake Ontario from the Canadian side, by 
extrapolating data from selected areas representative of 
larger hydrogeologic regions. This report, which describes 
the hydrogeology of the IFYGL Duffins Creek study area, is 
one in a series of several reports dealing with the ground^ 
water regimes of seven, selected representative areas 
along the Canadian shore of Lake Ontario. 



Toronto, May, 1979 G.H. Mills, Director 

Water Resources Branch 
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ABSTRACT 

The assessment of the ground-water contribution to 
Lake Ontario from the hydrogeologic region along the north 
shore of the lake, extending from the Humber River near the 
western limits of the City of Toronto to the City of Oshawa, 
was undertaken by the Ontario Ministry of the Environment 
as part of its contribution to the International Field Year 
for the Great Lakes (IFYGL) program. Data were extrapolated 
from the hydrogeologiG evaluation of the lower part of the 
Duff ins Creek drainage basin which was considered to be 
representative of the ground-water regiire in this hydrogeo- 
logic region. 

Field investigations of the geology, test drilling, 
data from water-well records on file with the Ministry 
and published literature were utilized for the assessment 
of the hydrogeology . 

Ground-water supplies in the study area are obtained 
from the upper part of the. Upper Ordovician shales and the 
overlying Wisconsinan drift. High-level, glacial lake 
deposits of permeable materials that are being utilized for 
grovmd water in the study area were tentatively identified 
as being part of the Thomcliffe Formation. Water-wells 
completed in the overburden are more productive than wells 
completed in the bedrock, with overburden wells having a 
specific-capacity range of 0.8 to 5.0 Igpm per foot of 
drawdown . 

An excess of 4,29 inches of precipitation over the 
long-term annual average during the IFYGL period (April 1, 
19 72 to March 31, 1973) resulted in an increase in storage 
of 0.85 feet in the ground-water reservoir of the study 
area. 

Permeability values, derived from short-term pumping 
tests, were assigned to the various materials in the study 
area. From these values it was estimated that five cfs of 
ground water is discharging directly into Lake Ontario 
along the shore of the hydrogeologic region. 
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Figure 1. Location of th« IFYGL Duffins Crattk study area in southern Ontario. 



INTRODUCTION 

PURPOSE & SCOPE 

Under the sponsorship of the Canadian and U.S. National 
cominittees for the International Hydrological Decade (IHD) 
program, the International Field Year for the Great Lakes 
program (IFYGL) was established. The purpose of the IFYGL 
was to study, in detail, the various hydrologic aspects 
associated with Lake Ontario and its drainage basin. Part 
of the Ontario Ministry of the Environment's contribution 
to the IFYGL program was the study of ground-water inflow to 
Lake Ontario, by extrapolating data from selected areas 
representative of larger hydrogeologic regions. 

The ground-water regime developed in any region is a 
result of the geology, topography, drainage and climate of 
that area. In the hydrogeologic region along the north shore 
of Lake Ontario, from the Humber River near the western limits 
of the dty of Toronto to the City of Oshawa (Figure 1.), 
ground-water conditions are considered to be uniform. The 
Duffins Creek drainage basin is situated in this area and 
was selected for the IFYGL study as a result of hydrogeological 
investigations being conducted by the Ministry in the basin. 
These investigations are part of an ongoing program for the 
assessment of water resources in the Province of Ontario. 

For the purposes of the IFYGL study, to determine the 
amount of ground-water inflow to Lake Ontario, only the area 
adjacent to the Lake was studied. Field investigations were 
made of the geology and test drilling was conducted to provide 
information on aquifer characteristics, subsurface geology 
and ground-water levels in the area. Other aspects of the 
hydrogeologic assessment utilized information obtained from 
published literature and water-well records filed with the 
Ministry of the Environment. 

LOCATION 

The study area is situated in the Regional Municipality of 
Durham, formerly the -Township of Pickerina in the County of Ontario, 
between longitudes 78 58' 15" W and 79^ 06' 00"W on the north shore 
of^Lake Ontario. At the north end, the area is bounded by latitude 
43 54" 30" N, Duffins Creek enters the north oart of the study 
area at longitude 79 04' 08" W and flows through the western part 
of the Town of Ajax into the Lake at longitude 79° 02' 08" W 
(Figure 1) , 
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GEOGRAPHY 



PHYSIOGRAPHY 

Two major physiographic regions, the Iroquois Plain 
and the South Slope as delineated by Chapman and Putnam (1966) | 
are found in the study area (Map 1) . All but the northwest 
comer of the study area is considered to be part of the 
Iroquois Plain. This plain has been defined by Chapman and Putnam 
as the lowland that borders Lake Ontario which was inundated in 
late Pleistocene times by glacial Lake Iroquois. In the study 
area, the lake bottom is ground moraine which was smoothed by 
wave action and the deposition of lacustrine sediments in 
depressions within the inundated area. The most prominent 
feature in the region is the northeast- trending bluff of the 
Iroquois shoreline (Map 1) , separating the Iroquois Plain in the 
south from the drumlinized till plain of the South Slope region 
to the northwest. 

A maximum relief of approximately 450 feet is present 
from Lake Ontario at + 245 feet above sea level (asl) to the 
highest spot elevation in the northwest corner of the study area 
at approximately 700 feet asl. A regional topographic gradient 
of approximately 25 feet per mile , extending inland from the 
Lake Ontario shore and increasing northwards, is present in the 
bulk of the area. This is the offshore facies of glacial 
Lake Iroquois where silt and clay deposits are found in the 
depressions of the inundated and wave-modified ground moraine,, 
In the western part of the study area, in the vicinity of 
Frenchman's Bay, the regional topographic gradient increases to 
75 feet per mile to the northwest. 

In tile nearshore facies of Lake Iroquois , the regional 
topographic gradient increases in the order of 100 feet per 
mile to the northwest. This area is delineated by boulder 
pavements and beach and bar deposits extending one to one and 
a half miles south of the Iroquois shoreline. In the drumlinized 
till plain of the South Slope physiographic region, north of 
the old shoreline, the regional topographic gradient increases 
to the northwest at approximately 150 feet per mile. 

Minor streams in the southwest comer and far eastern 
part of the study area drain directly into Lake Ontario. 
The remaining and major part of the area is drained by the lower 
reaches of Duffins Creek and its tributaries. The Duff ins Creek 
watershed is drained by two major creeks, the East and West 
branches. Their confluence is about three and one-half miles 
north of Lake Ontario near the northwest limits of the Town of 
Pickering, in the central part of the study area (Map 1) . 
The East Branch enters the study area in the north and flows 
southward, with an average fall of 30 feet per mile , for approxi- 
mately four miles to its junction with the West Branch. The 
latter enters the area from the west, flowing eastward for three 
miles, with an average droD of 40 feet per mile. A reduced 
gradient of less than 20 feet per mile occurs in the streambed 
from the union of the branches to the creek's mouth at Lake Ontario 
This reduced gradient, in the lower reaches of Duffins Creek, 
is in part a result of differential uplift causing partial 
drowning of the creek mouth. 



CLIMATE 

The climate of the study area is influenced by the 
presence of Lake Ontario. Most of the study area is in the 
climatic region of the "Lake Ontario Shore", as classified by 
the Ontario Department of Agriculture and Food (1966) , A 
small portion, in the north and northwest of the area, falls 
into the climatic region of the "South Slopes". 

Three precipitation stations. Frenchman's Bay, Pickering, 
Greenwood MTRCA and one meteorological station, Pickering- 
Audley , were in operation within the vicinity of the study 
area (Map la) . The meteorological data at the Pickering-Audley 
station were assumed to indicate the climate in the area. 
The long-term annual precipitation (11 years record) at this 
station is 31.91 inches and the normal long-term annual 
teitperature is 44.0°F. The normal monthly precipitation and 
temperature values for the Pickering-Audley station are shown 
in Table 1. These long-term data are not available for the 
Frenchman's Bay, Pickering and Greenwood MTRCA stations. 



TABLE 1. NORMAL MONTHLY PRECIPITATION AND TEMPERATURE DATA AT THE PIOCERING-AUDLEY 

METEOROLOGIC STATION (from Environment Canada, Atmospheric Environment Services, 1971). 
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19.6 


20.4 
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GEOLOGY 
BEDROCK GEOLOGY 

The study area is underlain by shales from the lower 
part of the Upper Ordovician time scale. These shales are 
part of the middle member of the Whitby Formation 
(Liberty, 1969) form.erly known as the Blue Mountain, 
Gloucester and Colling\^ood formations. The middle member 
is comprised of thin-bedded, brown, bituminous, micaceous, 
slightly calcareous shale which contains m.inor amounts of 
black, blue and dark carbonaceous shale. According to 
Liberty (1969) , these deposits are indicative of slow sedi- 
mentation in quite, shallow waters which is further supported 
by field evidence of ripple marks and abundant mud- cracked 
surfaces . 

The structure of the bedrock is obscured by a veneer 
of glacial drift in the study area. The bedrock is relatively 
flat-lying and dips to the south -southwest at approximately 
20 to 30 feet per mile (Liberty, 1969). 

The elevation of the bedrock surface is illustrated 
in Map 2. The attitude of this surface was compiled from 
records of water wells (on file with the Ministry of the 
Environment) that penetrate the overburden to the bedrock; 
borehole data obtained for the foundation investigation of 
the Pickering Generating Station (The Hydro-Electric Power 
Commission of Ontario, 1965) ; and foundation reports for 
provincial highways in the study area (Ontario Ministry of 
Transportation & Communications 1959, 1967, 1969 and 1970). 
The location of these data points are shown in Map 3. The 
data suggest that the bedrock surface rises gently northwards 
with an average gradient of 25 feet per mile, similar to that 
of the ground surface. The average elevation of the bedrock 
surface at the Lake Ontario shore is approximately 200 feet 
asl (50 feet below lake level) and 350 feet asl (25 to 100 
feet below the ground surface) in the north end of the study 
area. The maximum recorded relief is in the order of 200 
feet, approximately one half that of the topographic surface, 
from 14 5 feet asl near 'Lake Ontario at the Town of Ajax to 
350 feet asl at the north end of the study area. 

Two major bedrock valleys which drained to the south 
have been identified in the study area (Map 2) , The East 
Branch of the Duffins Creek occupies the same general area 
as the largest bedrock valley. Information to date suggests 
that tlie mouth of Duffins Creek at Lake Ontario is approxi- 
mately two miles west of this ancient valley. The other 
bedrock valley V7as occupied by an ancient stream that drained 
into the present-day Frenchman's Bay area from the northwest. 



PLEISTOCENE GEOLOGY 



Introduction 



The record of geological history from the Ordovician 
to the Quaternary periods is missing in the southcentral 
part of Ontario. The unconsolidated materials overlying the 
bedrock in the study area were deposited during the 
Pleistocene or Glacial Epoch and are generally considered to 
be Wis cons inan in age. Investigations of the glacial deposits 
have been made by earlier workers, notably Colem.an (19 36) , 
Caley et al (1947) , Olding et al (1956) , Gravenor (1959) , 
Chapman and Putnam (1966) and Hewitt (1969) . 

The deposits in the study area are described as glacial, 
glaciofluvial and glaciolacus trine in origin; with thick- 
nesses ranging from a few feet to more than 150 feet (Map 4) . 
In general, low areas in the bedrock surface, such as ancient 
stream valleys cut into the bedrock (Map 2) , have been filled 
with glacial debris; whereas bedrock highs (east and west of 
Duffins Creek near the shore of Lake Ontario) are covered 
with a relatively thin cap of glacial material. It would 
appear that the net result of glaciation in the study area 
was to smooth the pre-existing bedrock surface. A paucity 
of data on the bedrock surface in the northwest part of the 
study area does not allov? for accurate delineation of 
the thickness of the overburden deposits there. 

Stratigraphy and History 

Stratigraphy . • . . The most recent work on the 
Pleistocene stratigraphy, to the west and east of the study 
area, has been completed by Karrow (1967) and Singer (1973), 
respectively (Table 2) . For the most part, only the deposits 
of the Fourth (most recent) Pleistocene glaciation, the 
Wisconsinan Stage, are present in this part of Ontario. 
Some exposures of Illinoian Till from the Third Pleistocene 
glaciation and interglacial deposits of the Sangaraonian 
Stage have been identified in the Toronto area to the west 
of the study area. These deposits have not been identified 
east of Toronto* 

Early Wisconsinan and Older .•.. To date , depos i ts 
older than the Middle Wisconsinan have not been identified 
in surface exposures within the study area. Subsurface 
information, obtained from water-well records on file with 
the Ministry, suggests that the Scarborough Formation and 
Sunnybrook Till of Early Wisconsinan age may be present in 
the western part of the area. Positive correlation of 
these deposits has not been successful to the present time. 
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Scarborough Formation 
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- - 
Sangainonian 






Don Formation 




Illinoian 






York Till 




g 


Cincinna- 
ti an 


Edenian 


Gloucesteriar 


Nott^awasaga Group 




Whitty Formation 





The type section of the Scarborough Formation, as 
exposed in the Scarborough Bluffs, consists of a lower clay 
unit averaging 90 feet of exposure above the present lake 
level and an upper sand unit averaging 50 feet in thickness 
(Karrow, 196 7). Subsurface deposits west of Frenchman's Bay 
and east of the Town of Ajax may be correlative with the 
Scarborough Formation. 

In the Toronto area, the Sunnybrook Till has been 
established as a major till sheet extending under most of 
the area (Karrow, 196 7) . According to Karrow, the Sunnybrook 
Till is a fine-grained till varying from a silt till to a 
silty clay till with a low pebble content. The tiiickness 
of the till is variable and is most often found to be 20 to 
30 feet thick in the Scarborough Bluffs. Data from water- 
well records suggest that the Sunnybrook Till is present 
at depth throughout most of the study area. To the east, 
a basal till with similar characteristics has been identified 
by Singer (19 73) in the Bowmanville area and correlated 
with the Sunnybrook Till, 

According to Karrow (196 7) , the sediments of the 
Scarborough Formation are considered to be part of a delta 
formed in the Scarborough area at the mouth of a large 
southward-flowing river. This river eitptied into a 
proglacial lake (Lake Scarborough) that occupied tlie Lake 
Ontario basin when the earliest Wisconsinan glacier reached 
the eastern end of present-day Lake Ontario from the north. 
The ice dammed up the St. Lawrence valley raising the water 
level while drainage was diverted down the Mohawk and 
Hudson valleys to the sea. The ice front then retreated 
allowing the lake level to fall during which time the 
Scarborough delta was dissected by streams. A cooling of 
the climate brought about a major ice advance which 
completely covered the area and deposited the Sunnybrook 
Till. The ice at this time flowed into the Lake Ontario 
basin and spread westward over the Scarborough delta. The 
texture of this till is indicative of the fine-grained 
lacustrine material that was incorporated during the advance 
over the Lake Scarborough deposits in the Lake Ontario basin* 

Thorncliffe Formation .. . Karrow (19 6 7) has named 
the variety of glacial, fluvial and lacustrine deposits 
found between the earliest Wisconsinan till, the Sunnybrook, 
and the latest till, the Leaside, the Thorncliffe Formation, 
In Karrow' s type section on the Scarborough Bluffs, over one- 
half the sediment is stratified fine sand, the rest being 
dense silt and varved clay. Karrow indicates that the 
genesis of these deposits appears to be lacustrine with 
stream and deltaic sediments occurring at higher elevations. 



The Thorncliffe Formation has been identified by 
Karrow (196 7) in exposures along the West Branch of Duff ins 
Creek. For the present study, these exposures have been 
correlated with subsurface data available from water-well 
records (Map 3) on file with the Ministry of the Environment. 
The linear distribution of permeable materials (sands and 
gravels) in the sxjbsurface and their relative position to 
one another, as outlined on Map 5, suggests that a series 
of encroaching, high-level lakes were initiated after the 
withdrawal of the Early Wisconsinan ice sheet that deposited 
the Sunnybrook Till. These high-level proglacial lakes 
were formed by a readvance of the ice into the Ontario basin, 
with a probable ice front near the middle of the basin. 
Three major phases of lake levels, early, intermediate and 
late, have been identified from the distribution of 
sediments in the subsurface at elevations of 245 to 301 
feet above sea level (asl) , 298 to 347 feet asl and 348 to 
425 asl, respectively (Maps 6 to 9 inclusive) . Associated 
features such as bars, spits and deltas appear to have 
existed during the various lake stages. Iliese beach 
deposits have been interpreted to be part of the Thorncliffe 
Formation in the study area, although their complex 
nature is only poorly landerstood due to the lack of control. 
Some similar, deltaic- like deposits (the Clarke Deposits) 
have been identified at an elevation of approximately 320 
feet asl in the Lake Ontario bluffs near Bowmanville by 
Singer (19 7 3) . It is suggested that the Clarke Deposits 
in the Bowmanville area are equivalent to the early lake 
phase of the niomcliffe Formation at an elevation of 245 
to 301 feet asl in the Duffins Creek area; assuming that the 
westward decrease in elevation due to differential uplift 
is approximately 30 feet as indicated by the tilt in the 
Lake Iroquois shoreline* 

Seminairy Till . . . According to Karrow (196 7) the 
Seminary Till is a silty sand till varying in thickness 
from nearly zero up to 26 feet and represents a minor 
advance of no great inland extension in the Scarborough 
area. Surface exposures of the Seminary Till have not 
been identified in the study area. Subsurface information 
suggests that this till is present and appears to represent 
an ice advance between the early and intermediate lake 
phases of the Thorncliffe Formation. 

The successive shorelines of the intermediate lake 
phase of the thorncliffe Formation, as represented by the 
deposits that overlie the Seminary Till, generally conform 
to that of the early lake phase but at a higher elevation. 
This suggests that the ice advance that deposited the 
Seminary Till did not modify the existing landscape 
to any great extent and did not retreat to the previous 
ice-front position. The sandy nature of this till may 
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be in part attributed to overriding and incorporation of 
the beach sand from the early lake-phase deposits . 

Meadowcliffe Till . . . . The Meadowcliffe Till has been 
tentatively identified at the bottom of an abandoned gravel 
pit near the western limits of the Town of Ajax^ east of 
Duff ins Creek. The appearance of this till is similar to that 
of the Sunnybrook Till although carbonate analyses of the 
matrix material show differences as illustrated in Table 3 
below I 

TABLE 3. COMPARISON OF THE MEADOWCLIFFE AND SUNNYBROOK 
TILLS (Harrow, 1967) 



Till Total Carbonate (%) Cal cite/Dolomite 



Sunnybrook Till 12% Qti 

(22 samples) 

Meadowcliffe Till M. 4i.i 

(9 samples) 



^toalysis of the till from the gravel pit in the Town of Ajax 
indicated a total carbonate content of 32.5% and a calcite to 
dolomite ratio of 4.81. 

In the Scarborough Bluffs, the Meadowcliffe Till has a 
uniform thickness of 35 to 40 feet and reaches up to a thick- 
ness of 80 feet inland. Karrow (1967) indicates that this 
till appears to extend farther to the north than the under- 
lying Seminary Till, but is still of limited, known extent. 
Singer (19 7 3) has identified the Meadowcliffe Till in the 
Lake Ontario bluffs in the Bowman ville area. Only a minimum 
of information is available regarding the Meadowcliffe Till 
in the study area and data from water-well records suggest 
that this till likely represents an ice advance between the 
intermediate and late, lake phases of the Thorncliffe 
Formation, 

Karrow (1967) indicates that the fine texture of the 
Meadowcliffe Till is a result of the incorporation of fine- 
grained lake sediments that were deposited in a proglacial 
lake that was created when the ice sheet that deposited the 
previous till sheet (Seminary Till) retreated into the 
Lake Ontario basin. Both the Seminary and Meadowcliffe 
Tills are considered to represent minor ice advances, 
indicating a general cooling of the climate just prior to 
the major Late Wisconsinan ice advance. 
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Leaside Till . . . The Leas ids Till represents the 
last Wis eons inan ice advance in the Toronto-Bowmanville 
area. Where exposed (Map 4) , this till is observed to be 

a silty sand till with a high stone content. Carbonate 
analyses of the till matrix in the study area are com-- 
parable to those values reported by Karrow (19 67) in the 
Scarborough area and are listed in Table 4« 



TaBLE 4, CARBONATE ANALYSES OF THE LEASIDE TILL 



^rea Total Carbonate (%) Calcite/Dolomite 

Scarborough (Karrow> 1967} 31.0 4.0 

(26 samples) 

Duffins Creek 30.6 |.ai 

(9 samples) 



A break in the deposition of the Leaside Till^ 
represented by a thin sand seam or discontinuous beds of 
sand or sand and gravel up to 18 feet in thickness » has been 
observed in some sections along the present-day Lake Ontario 
bluffs and along the valley walls of Duffins Creek. AcGording 
to Karrow (196 7) , these two till sheets represent two ice 
advances separated by a very short time interval and are 
considered to be part of the Leaside Till. This same feature 
was observed in the Bowroanville area by Singer (1973) . 

In some sections along the Lake Ontario bluffs , up to « 
seven feet of the uppermost Leaside Till shows (one-half inch 
thick) thin pseudo -bedding,, grading upward from a silty sand 
till, through a silty clay till with few stones into beds of 
varved clay. This feature is probably indicative of the 
retreat of the ice edge in standing water that was dammed up 
around the periphery of the ice. 

Lake Iroquois Deposits^ . . . Karrow (196 7) indicates 
that high-level lakes were formed during the retreat of the 
ice front to the eastern end of Lake Ontario, Meltwaters 
were dammed up between the higher land surface and the ice 
front that was retreating to the east in the L^e Ontario 
basin. The levels of these peripheral lakes dropped as 
lower outlets were uncovered by tiie retreating ice. Lake 
Iroquois was formed when an outlet stabilized at Rome, 
;New York, while the ice retreated to the St, Lawrence valley. 



Wave action in Lake Iroquois modified the drumlinized 
till plain that resulted from the ice sheet that had 
deposited the Leaside Till over the area. A wave -eroded 
bluff, (Map 4) ranging in height from 20 to 70 feet and 
trending northeast across the northwest comer of the study 
area, was cut into the till plain. Ihis prominent feature 
serves as a physiographic boundary separating the wave- 
modified till plain or lake plain to the southeast from the 
irregularly rolling and ice-fluted surface to the northwest. 
Longshore currents in Lake Iroquois removed the material in 
the bluffs and exposed older beds where the till cover was 
thin. The base of the wave-cut bluff is at an elevation 
of approximately 4 80 feet asl and is coincident in part 
with the shorelines from the late lake phase of the 
Thorncliffe Formation (Map 5) commencing at an elevation of 385 
feet above sea level. The writer believes that some of the 
sediments previously mapped as Lake Iroquois deposits belong 
to the Thorncliffe Formation and the surface materials may 
have only been reworked by currents in Lake Iroquois. 

Granular deposits of sands or sands and gravels 
constitute the near-shore fades of Lake Iroquois. A 
relatively- thin veneer of beach sands is generally found in 
close proximity to the Iroquois bluff. In some areas, where 
erosion was dominant, a boulder pavement derived from the 
winnowing of fines from the coarser material comprising the 
Leaside Till, marks the high -energy zone of Lake Iroquois. 
Within the study area, two embayments are present in the 
Iroquois shoreline across which gravel bars were built. 
These embayments were low areas between drumlin ridges that 
were carved by the glacier on the Leaside Till plain. Erosion 
of the headlands resulted in the development of high cliffs 
on the ends of the drumlins that projected into Lake Iroquois . 
The sediment from the cliff erosion was deposited as gravel 
bars across the bays. The present-day East Branch of Duffins 
Creek and a tributary to the west of the Creek breach the 
shoreline at these embayments (Map 4) . In contrast, the 
West Branch of Duffins Creek cuts through a straight stretch 
of shoreline, in a gorge approximately 150 feet deep, having eroded 
over 100 feet of shoreline since the disapoearance of 
Lake Iroquois. Much of the granular deposits have been 
exploited for building materials and only remnants of these 
deposits are present in the area today. 

Deep-water sedimentation partially filled the 
depressions in the offshore area covered by Lake Iroquois. 
Fine-grained deposits of varved clay and stratified silt 
tend to blanket the underlying till plain to the edge of 
the present-day Lake Ontario shore. Drumlin ridges, trending 
northwest-southeast in the southern part of the area, rise 
through the lake plain and do not appear to have been greatly 
modified by the wave action in Lake Iroquois; whereas some 
of the higher land to the north appears to have been 
levelled and smoothed. 



When the St. Lawrence valley became free of ice, 
the waters of Lake Iroquois were drained. Remnants of 
shorelines at 350 and 300 feet asl (Map 4) suggest that 
the lake level dropped in stages as the ice front retreated 
and lower outlets to the sea were uncovered. None of the 
shorelines are well-developed, suggesting that the time 
interval between the stages was of short duration. 
Karrow (1967) indicates that Lake Iroquois drained to more 
than 200 feet below the present level of Lake Ontario 
approximately 10,000 years ago. Differential uplift of 
the land surface has gradually raised the water in the 
Ontario basin to its present level at approximately 24 5 
feet asl and the formerly horizontal shoreline of 
Lake Iroquois has been tilted up to the northeast 
Coleman (19 37) reported that the shoreline was at an 
elevation of 362 feet asl in the City of Hamilton at the 
west end of Lake Ontario and 730 feet asl towards the 
eastern end of the Lake , at Pancake Hill north of the City 
of Belleville. The differential rise in the shoreline 
increases to the nortiieast along the direction of maximum . 
tilt at N 20° E (Coleman 1937). Karrow (1967) indicates 
that the average rise of the shoreline in the Scarborough 
area, adjoining the study area on the west, is approximately 
2.5 feet per mile. 

Recent Deposits . . . Alluvial deposits are found along 
most of the stream courses in the study area. The material 
forming the banks and underlying the streams is subjected to 
erosion whereby the finer material is carried downstream 
and deposited along the flood plain and streambed- The 
.coarser material, such as boulders in the till, is left 
behind to form a lag deposit in the streambed 

The lower reaches of the streams flowing into the 
lake are presently being submerged by the rising water 
level in Lake Ontario and marshes have developed near their 
mouths. Notable examples in the study area are the marshes 
around Frenchman's Bay and the mouth of Duff ins Creek (Map 4), 

The sediments along the Lake Ontario beach are the 
result of the erosion of the materials in the present-day 
bluffs and consist mainly of sand and some gravel blanketing 
older deposits. Large boulders, eroded out of the till, 
litter the beach in some areas and help to protect the base 
of the bluffs from erosion. 
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HYDROGEOLOGY 

GROUND-WMER OCCURRENCE 

Ground water is an aspect of the hydrologic cycle which 
has resulted from the infiltration of moisture at the earth's 
surface, moving through the "unsaturated zone" into the 
"saturated zone" . The water stored in the saturated zone is 
called ground water. Under the influence of gravity, ground 
water moves to a point of natural or artificial discharge in a stream, 
lake or well. Natural recharge to this system occurs either directly 
from precipitation infiltrating into the ground or indirectly 
from snowmelt , standing bodies of water and surface streams. 

An aquifer is defined as a water-bearing horizon that 
will transmit significant amounts of water into a well. In 
the exploration for ground-water resources, the connotation 
of the wording "significant amounts" further defines the term 
aquifer with respect to the usage of water for domestic, 
industrial or municipal supplies. An aquifer in the context 
of this report is one that will supply water for at least 
domestic use. In the study area, aquifers are found both 
within the overburden and in the bedrock. 

Ground water is stored in the interstices or openings 
and pore spaces between the grains comprising the material 
through which the ground water moves. The capacity of any 
material to absorb, yield or hold water is dependant on the 
porosity and permeability of the material (i.e. the total 
number of interstices or voids and the interconnection of these 
interstices, respectively). Thus, a high porosity material 
such as shale or clay will not yield large quantities of water 
due to absorption phenomena resulting from the molecular size 
and poor interconnection of the pore spaces. 

Within the overburden, the most productive materials 
are the granular deposits of sand and gravel. The extent 
of these granular deposits has been delineated in Maps 6 to 
9 inclusive, for the study area, as indicated from data in 
water-well records on file with the Ministry of the Environment. 
In a previous section of this report, it was suggested that 
these granular materials represent ancient beach deposits 
that were overridden by the latest Wisconsinan glaciation. 
These deposits are presently being utilized as the primary 
source of ground water in the study area. 

Numerous bored and dug wells penetrating low water-yielding 
materials such as glacial deposits of till or glaciolacustrine 
clays and silts are also being used as a source of ground 
water. The yield from the wells completed in these materials 
is essentially controlled by the storage capacity of the 
large-diameter wells. Recent studies in the prairies 
(Meyboom et al, 19 66) suggest that the movement and storage 
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of ground water in glaqial tills is within joints developed 
in the till. 

The shale bedrock in the study area is generally a poor 
aquifer. Microscopic examination indicates that its porosity 
is low resulting from the consolidation and compaction of 
the shale since its deposition. Joints, fractures, cleavage 
planes and faults were siJbsequently developed in the bedrock, 
probably as a result of tectonic processes. Other processes, 
such as those associated with weathering, have enlarged these 
planes (including the original bedding planes) creating 
secondary porosities through which the ground water presently 
flows. These secondary openings tend to become tighter and 
fewer in number with increasing depth. 

GROUND-WATER LEVELS 



The change in storage of the ground-water reservoir can 
be identified by the measurement of rising and declining water 
levels in water wells or boreholes. During periods of high 
precipitation, storage in the ground-water reservoir is 
increased, after the soil moisture deficiency (SMD) in the un- 
saturated zone is replenished. When the SMD is satisfied, the 
infiltrating water will percolate downwards into the saturated 
zone under the influence of gravity. Conversely, storage is 
decreased under conditions where the infiltrating water is 
intercepted before reaching the saturated zone, while discharge 
from that zone continues. These conditions occur during 
periods of low precipitation and/or excessive use of ground 
water by man or vegetation. 

An established, observation-well network of private and 
abandoned water wells (Sibul et al, 19 77) was supplemented 
with three observation-well installations constructed by the 
Ministry of the Environment for the IFYGL program (Map la) . 
The change in water levels was computed for the wells in this 
network over the Field Year period from April 1st, 19 72 to 
March 31st, 1973 and the results are presented in Table 5. 
The data in Table 5 suggest than an overall increase in the 
storage of the ground-water reservoir occurred during the 
Field Year. This is indicated by the average rise of 0.85 
feet in the water levels of the observation-well network from 
April 1, 1972 to March 31, 1973. 



TABLE 5. STATIC LEVELS AMD CHANGE IN WATER LEVELS 

(s) OF THE OBSERVATION -WELL NETIVORK IN THE 
IFYGL DUFFINS CREEK STUDY AREA 









Static 


Static 












Level 


Level 




Average 


Well 


Diameter 


Depth 
12' 


Apr. 1/72 
7.80' 


Mar. 31/73 
8.60' 


s 
-0.80 


s 


329 




308 


i 


14' 


3.34' 


3.35' 


-0.01 




P-2 


t«' 


16' 


6.19' 


5.85' 


+0.34 




P-IA 


2" 


17' 


4.88' 


2.75' 


+2.13 




P-IB 


2" 


17' 


5.23' 


2.87' 


+2.36 




302 


48" 


20 i 


2.90' 


0.60' 


+2 . 30 




301 


48" 


23' 


17.50' 


15.45'E 


+2.15 


+1.21 


P-IA 


h" 


30' 


4.64' 


2.84' 


+ 1.80 




P-2 


¥' 


36' 


7.80' 


5.98' 


+1.82 




330 


'^W. 


40' 


25.00 'E 


25.09' 


-0.09 




303 


48" 


44' 


5.04' 


4.58' 


+0.46 




P-IA 


h" 


45' 


4.66' 


3.17' 


+1.49 




336 


^w 


47' 


40.80'E 


40.17' 


+0.6 3 




P-2 


^'^ 


49' 


6.53' 


7.27' 


-0.74 


+0.62 


331 


2" 


73' 


28.70'E 


28.53' 


+0 . 17 




337 


2" 


92' 


39.92'E 


38.26' 


+1.65 




332 


2" 


157' 


13.80'E 


13.62' 


+0.18 


.- 


333 


2" 


238' 


12.95'E 


12.65' 


+0.30 




304 


6" 


424' 


189.43' 


189.38' 


+0.05 


+0.47 



+0 . 85 



E = estimated from observation-well hydrograph 



Climatological data, presented in Table 6 below, 
indicate that excess precipitation over the long-term annual 
average contributed to the rise of ground-water levels that 
occurred in the study area during the Field Year. It is assumed 
that the excess of 4.29 inches over the long-term annual 
average of 31.91 inches recorded at the Pickering-Audley 
station, occurred over all of the study area. This excess 
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precipitation is reflected by the overall rise of 0.85 feet 
in the ground-water levels of the observation-well network. 



TABLE 6 , 



AVERAGE ANNUAL PRECIPITATION IN THE IFYGL DUFFINS 
CREEK STUDY AREA 



Annual Preeip. 
Station (1/4/72-31/3/73) 

Frenchman's Bay 36.77" 

Pickering 35.82" 

Greenwood MTRC 35.58" 

Pickering- Audley 36 , 20 " 

Station Average 36. Of" 



Ave r a ge Lon g ~Te rm 
(12 yrs) Ann.Preclp. 



31.91" 



Deviation 
from Long 
Term 



+4.29" 



The water-level data (Table 5) suggest that the 
magnitude of the overall rise in water levels for the Field Year 
period decreases with increasing depth of the well. Wells 
under 30 feet in depth show an average increase in water 
level of 1.21 feet, whereas wells beb^een depths of 30 to 50 
feet show a smaller increase of 0,62 feet. This feature is 
also illustrated in site Pl-A near the shore of Lake Ontario, 
where observation-well installations at depths of 17, 30 and 
45 feet showed a progressive decrease in the net water-level 
change of +2.13, +1,80 and +1.49 feet, respectively, over the 
same time period. 

In addition, it appears that water levels at site 
Pl-A are influenced by changes in the level of tlie Lake as 
indicated by the similarity of the lake and well hydrographs 
for the Field Year period (Figure 2) . The high and low -water 
stages for both Lake Ontario and the observation wells at 
site Pl-A occurred during the months of June and November, 
respectively (Figure 2) , whereas the high and low-water staaes 
of the other observation wells (inland from Lake Ontario) 
occurred during the months of April and September, respectively 
(Figure 3) . 

GROUND -WATER MOVEMENT 
Ground-Water Flow Field 



The ground-water flow field is a three-dimensional 
field that occurs from the top of the saturated zone down to 
a depth where ground-water flow becomes impossible (i.e. where 
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FtgMfe 2. Hydrograph comparison of Lake Ontario and observation wells at site PI -A in the IFYGL Duff ins Creek studv area. 
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Observation welt 301 

Total depth: 23 feet ( 7 .01 m) 
Surface elevation: 735 ft amsl 
Well diameter: 48 inches (121 92cm) 
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Observation well 302 

Total depth: 20 feet (6.096m) 
Sur[3ce elevation: 830 ft. amsl 
Well diameter: 48 inches (121.92cm) 
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Observation well 303 

Total depth: 44 feet (l 3 . 4 1 1 m) 

Surface elevation: 325 ft. amsl 

Well diameter: 48 inches (l 21.92cm) 
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Figure 3. Observation-well hydrograplis for wells 301. 302 and 303 in the vicinity of the IFYGL Ouffins Creek study area. 



the material in the saturated zone becomes completely 
impenrieable) . The ground-water head at any point in this 
field can be expressed in terms of total head defined by 
Hubbert (19 5 3) as the siim of pressure head, elevation head 
and velocity head at that point; however, as the velocity 
head is so small for ground-water considerations it can be 
neglected. The static water level (the level to which water 
rises) in a non-pumping well is a measure of the groiand-water 
head at the point where the water was found. 

The total- flow vector (Meyboom et al, 19 66) of the 
three-dimensional ground-water field can be resolved into 
two, two-dimensional horizontal and vertical vectors. From 
these, the "flow resultant" is the flow component in the 
vertical plane and the "horizontal component" is the flow 
component in the horizontal plane (Meyboom et al, 19 66) , In 
either case, the component of flow that is depicted ("flow 
resultant" or horizontal component") is assumed to essentially 
represent the total-flow vector of the ground-water flow field. 

Water-Level Data 

As indicated above, the static water levels in wells 
are a measure of the potential 'causing ground-water flow. 
Water-level data from water-well records on file with the 
Ministry of the Environment indicate that the average static 
water level for 342 overburden wells in the study area is 15.5 
feet below the ground surface with a standard deviation of 
12.6 feet. These data suggest that the static water levels of 
wells completed in different water-bearing horizons within the 
overburden are usually found within a few feet of the land 
surface, regardless of the season or year in which the wells 
were drilled. Additional information from the observation- 
well network (Sibul et al, 19 77) indicates that the maximum 
fluctuation of ground-water levels has been less than 10 feet 
over the past three years of record (Figure 3) , or less than 
one standard deviation from the mean, static water level 
indicated from the water-well records on file with the Ministry 
of the Environment. The similarity in the data suggest that the 
general direction of ground-water flow can be deterrined from 
these water-level data. Anomalously deep water levels can usually 
be traced to excessive withdrawals or consumption of ground 
water or to deep wells completed in the less active or stagnant 
part of tiie ground-water reservoir, 

Potentiometric Surface 

The water-level data from water-well records in the 
study area were used to compile a potent iometric-sur face map 
for the ground water in the overburden (Map 10) . A potentio- 
metric map was not prepared for wells completed in the bedrock 
due to the paucity of wells penetrating the bedrock in the 
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studly area. The hydraulic-head values (i.e. water-level 
elevations) for these wells are shown separately on the map 
prepared for the overburclen (Map 10) . The similarity in 
head values for the two sets of data, overburden and bedrock, 
suggests that they are hydraulically connected. For the 
purposes of this report, the overburden and uppermost part 
of the bedrock are considered to constitute the active part 
of the ground-water reservoir. 

With respect to the potentiometric-surface map 
{Map 10), lines of equal potential were interpolated between 
values of hydraulic head obtained from the vjater-well records 
on file with the Ministry of the Environment. By definition, 
ground-water flow is at a direction normal to the lines of 
equal potential. Examination of Map 10 indicates that the 
potentiometric surface conforms to the topographic expression 
of tlie ground surface and may be considered as a subdued 
replica of the physiography in the study area. The ground- 
water divides appear to be generally coincident with the 
surface-vjater divides and by analogy, areas of recharge and 
discharge are topographic highs and low, respectively . Thus 
ground-water flow, in a lateral direction, is controlled by 
the topography; movem.ent being from topographically high to 
topographically low areas . 

Hydraulic Gradient 

Assuming that the "horizontal component" of flow 
(as illustrated in Map 10) represents the total-flow vector 
of the ground-water flow field, the slope of the potentiometric 
surface is the hydraulic gradient under which ground-water 
movement occurs. Ground-water movement is dependant on the 
hydraulic gradient according to Darcy's Law, Q = KIA , where: 

Q = quantity of flow 

K = coefficient of permeability of the porous material 
I = hydraulic gradient 

A = cross -sectional area through which movement takes 
place . 

The hydraulic gradient (I) is the head loss (in feet of water) 
per unit distance (foot) of flow path, resulting from the 
frictional resistance to flow within the voids of the material 
through which the ground water moves. As indicated earlier, 
the potentiometric surface (Map 10) is similar to the land 
surface and ground-water flow is from high to low potential 
in a direction normal to the lines of equal potential. The 
hydraulic gradient (I) is in the order of 0,005 (feet per foot) 
for the bulk of the study area to 0.014 (feet per foot) in the 
vicinity of Frenchman's Bay, near the western part of the study 
area . 
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If the "horizontal component" of ground-water flow 
does not represent the total-flow vector of the ground-water 
flow field, then the slope of the po ten tiome trie surface is 
only a component of the true hydraulic gradient. It would 
then be assumed that the "flow resultant" (in the vertical 
plane) represents the total-flow vector of the ground-water 
flow field. In areas where there is a strong downward or 
upward component of flow, the vertical gradient determined 
from observation wells or piezometers, can approximately 
represent the hydraulic gradient. The vertical gradient is 
the head loss (in feet of water) per unit distance (foot) of 
depth. 

For the Field Year program, piezometer nests were 
installed at several sites (Map la) in order to determine the 

potential distribution of the ground water in the subsurface * 
Head measurements (static water levels) , obtained from these 
piezometer installations at depth, were used to compute the 
vertical gradients at the sites and the results are presented 
in Table 7 below: 



TABLE 7. COMPUTED^ VERTICAL HEAD GRADIENTS IN PIEZOMETER 
NETS Pl-A, Pl-B and P2 ON MAY 2, 19 72, 
NOVEMBER 6, 1972 and MARCH 29, 1973. 



Nest Material De 



PI- A overburden 17' 

overburden 30' ' 

overburden 45' 



Pl-l overburden %7 ' 



overburden 31' 



P2 overburden 16' 

overburden 36' 

bedrock 49 ' 

*NA = data not available. 



-0.020 



Vertical Gradient 



2/5/72 6/11/72 29/3/7 3 Average: 

+0.00 7 +0.020 ^0.00 7 +0.00 7 
-0.00 7 +0.016 ^0.022 -0.004 



lA* 



NA* -0 . 02 



-0.011 -0,021 -0.007 -0.013 



"0.101 -0.117 -O.IOQ -0.12 3 



m^ 



The vertical gradient ranges from -0.167 to +0.020 
(feet per foot of depth) in the study area, as illustrated 
in Table 7 above. Intermittent upward movement is indicated 
at site Pl-A, near Lake Ontario, at various times of the year, 
suggesting that this is an area of ground-water discharge. 
The negative vertical gradients measured at other times of 
the year may be local effects due to fluctuating water levels 
in Lake Ontario. Average vertical gradients at site Pl-A 
have the same order of magnitude as the estimated hydraulic 
gradient for the "horizontal component" of ground-water flow 
in Map 10. This suggests that the total flow-vector of the 
ground-water flow field can be approximated by either the 
vertical or the horizontal vector of flow at site Pl-A (ie 
"flow resultant" or "horizontal component"). 

The negative vertical gradients computed for site P2 
(Table 7) suggest that downward moveii^nt of ground water occ- 
urred at this site throughout the period of investigation. 
Within the overburden at this site, the vertical gradient is 
similar to the estimated hydraulic gradient for the horizontal 
component of ground-water flow (Map 10). This suggests again, 
that either the "horizontal component" or the "flow resultant" 
approximates the total-flow vector of the ground-water flow 
field. Strong downward movement of ground water is indicated 
by the average vertical gradient of -0.12 3 between the over- 
burden and bedrock at this site. If natural ground-water flow 
is wholly in the vertical direction at this site, the observed 
vertical gradient, between the overburden and the bedrock, 
must reflect a decrease in permeability of the bedrock (i.e. 
under these conditions the product of gradient and permeability 
must remain constant from Darcy's Law, Q = KIA , indicating that 
the permeability of the bedrock is less than that of the over- 
burden) . 

In summary, the hydraulic gradient under which ground- 
water movement occurs is similar to the regional topographic 
gradient which ranges from 0.005 to 0.014 feet per foot. The 
area near the Lake Ontario shore appears to be in a zone of 
general ground-water discharge. The area inland from the 
Lake Ontario shore appears to be in a general area of recharge, 
including those areas close to creek valleys (e.g. site P2) . 

AQUIFER CHAEACTERISTICS 

Introduction 



The water-yielding properties of an aquifer are a function 
of pore size and their interconnection. For quantitative^ 
ground-water studies, the most common method for determining 
aquifer characteristics is by means of a pumping test. The yield 
of a well, as determined from a pumping test, may be used to 
provide information on the water-yielding properties of the 
various aquifers encountered. In addition, pumping- test data may 
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be used to provide some estimate of the coefficients of 

trans mis sibility and permeability ( common ly referred to as ^ 

transmissibility and permeability/ respectively) of the 

materials through which the ground water is flowing. The 

coefficient of transmissibility (T) is defined as the rate 

of flow of water in Imperial gallons per day (Igpd) through 

a vertical strip of the aquifer one foot in width and 

extending the full saturated thickness of the aquifer under 

a hydraulic gradient of one foot per foot at the prevailing * 

temperature of water. The coefficient of permeability (K) 

may be determined using the relation -' 

jj _ T, where : 

m 

K = hydraulic conductivity (coefficient of per- 
meability, Igpd/ ft ) . 

T - coefficient of transmissibility (Igpd/ft) . 

m = contributing thickness of the aquifer (feet) . 

In order to evaluate the quantity of ground water discharging 
to Lake Ontario, an estimate of the permeability or trans- 
missibility of the materials through which the ground water I 
flows was made using well-yield data. 

Well Yields 

Statistical analyses of specific-capacity data (derived 
from pumping tests) by means of specific-capacity frequency 
graphs (Figure 4) provide a convenient method of comparing 
the water-yielding properties of the different aquifers in 
the study area. The specific capacity of a well is its yield, 
in Imperial gallons per minute per foot of drawdown (Igpm/ft) , 
for a stated pumping period and rate , and is numerically expressed 
as Q/s , where : 

Q = pumping rate in Imperial gallons per minute (Igpm) 
s = drawdown (loss of head) in feet. 

The specific-capacity data for wells completed in the 
overburden and bedrock were tabulated in order of magnitude .-| 

and the frequencies were computed by the Kimball (1946) method, | 

Values of specific capacity were plotted against percent of 1 

wells on logarithmic-probability paper. The slope of the line | 

indicates the variability of the water-yielding properties of | 

the material; steeper slopes signify greater variability than J 

flatter slopes. Comparison of the frequency curves for the | 

overburden and the bedrock (Figure 4) indicate a similarity in | 

slope (i.e. variability), with overburden wells being more 
productive than wells completed in the bedrock. For example, 
50 percent of the wells completed in the overburden have a 
specific capacity of 0.5 Igpm per foot or more, whereas 50 
percent of the wells completed in the bedrock have a specific 
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Figure 4, Specific-capacity frequency graphs for bedrock and overburden wells in the 
IFYGL Duff ins Creek study area. 
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capacity of 0.1 Igpni per foot or more. The specific 
capacities of overburden wells generally range from 0.0 8 to 
5.0 Igpm per foot of drawdown and the specific capaGities 

of bedrock wells range from 0.0 7 to 0.7 Igpm per foot of 
drawdown. 

The overburden freqiiency curve eons is te of three dis- 
tintive segments of different slopes. Two of these changes 
in slope appear to be restricted to the higher-yield gravel 
wells, suggesting that gravel wells with a specific capacity 
of 1 to 10 are more variable than those wells above or below 
that range. A change in slope also occurs in the range of 0.4 to 
0.7 Igpm per foot of drawdown for the bedrock frequency curve 
with slopes on either side of this range being the same. 
"Hie reason for this slope variability has not been determined 
to date; however, it may be a result of increased contribution from 
the overburden to the well through faulty well construction 
or artificial fracturing of the bedrock surface during drilling, 
permitting some leakage of the water from the overburden into 
ttie uncased portion of the well in bedrock, 

Transmissibility (modified non-equilibrium well formula) 

Information from water-well records on file with the Ministry 
of the Environment were utilized to estimate the transmissi- 
bility (and by analogy the permeability) of the different 
materials from short-term pumping tests of wells in the study 
area. Transmissibilities were obtained using the modified 
non-equilibrium well formula (Cooper and Jacob, 19 64), 

fp , 264 Q , where: 
s 

T = coefficient of transmissibility (Igpd/ft.) 

Q = pumping rate in Imperial gallons per minute (Igpm) 

s = change in drawdown per log,Q cycle (feet) 

The derivation of the formula is based on the following 
assumptions (Johnson, 19 66) : 

II The water-bearing formation is uniform in character 
and permeability in both horizontal and vertical 
directions . 

2) The formation has uniform thickness. 

3) The formation has infinite areal extent, 

4) The formation receives no recharge from any source. 

5) The pumped well penetrates and receives water from the 
full thickness of the water-bearing formation. 

1} The water removed from storage is discharged instan- 
taneously with lowering of the head. 
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Unfortunately, due to the nature of the data, (lack of 
observation wells, partial penetration of the aquifer by 
wells and non-uniformity of the aquifer, etc.) , the above 
assumptions are not met in totality. It is assumed that the 
transmissibility values obtained from the well-record data are 
within an order of magnitude of the absolute values. 

The results of the computations of permeabilities 
for the various material in the study area are presented 
below (Table 8) . Permeabilities were conputed from trans- 
missibility divided by contributing aquifer thickness (T/m) . 

TABLE 8. PERMEABILITIES (k) OF THE MATERIALS IN THE IFYGL 

DUFFINS CKEEK STUDY AKEA 

Material Number of Mean K- Std. Dey . Range of 
Wells (gpd/ft^) (gpd/ft^) (gpd/ft^i 

Gravel 33 468 822 7 - 3950 

Sand & gravel 11 366 455 3 - 1480 

Sand 14 77 75 2 - 250 

Till, clay, silt 35 15 14 1 ~ 59 

Bedrock 42 fi 9 t*! ^ 35 



Transmissibility (specific capacity) 

The coefficient of transmissibility of an aquifer can be 
estimated from the specific capacity of production wells. The 
tiieoretical specific capacity of a well may be written as, 



s — 


T 


, where 1 




114-^ ^(u,rw/3) 


u = 


2242 rw^S 
Tt 




rw _ 


rw 





^ T/P'm' 

W/ /n\ = well function for leaky artesian aquifers 
(u,x^./Bj (Han tush iy56) 

Q = discharge of pumped well (Igpm) 
s = drawdown (feet) 
rw = nominal radius of well (feet) 



T = coeffleient of transmisaibility (Igpd/ft.) 

S - coefficient of storage 

t = ptamping period (minutes) 

p' = coefficient of vertical permeability of confining 
bed (Igpd/ft ) 

m' = saturated thickness of confining bed (feet) 

Reliable data are not available on the coefficient of storage 
(S) in the study area. The magnitude of S depends on the 
elasticity of the aquifer material and the fluid. In confined 
aquifers, S does not show large variations and ranges from 
10-6 to 10"^ (Kruseman et al 19 70) . Almost all the wells in 
the study area on file with the Ministry of the Environment 
indicate confined conditions, having a static water level 
above that where the water was found. According to Meyer 
(Bentall, 1963) any changes in the coefficient of storage (S) 
correspond to only small changes in the transmissibility (T) 
and specific capacity (Q/s) ; therefore, inaccuracy in estimating 
S generally is not a serious limiting jEactor. 

With regard to the vertical permeability of the confining 
bed (P ' ) , data have not been compiled on this parameter in 
the study area to date. 

In addition to being dependant on the hydraulic proper- 
ties of the aquifer and confining beds, the specific capacity 
of a well will vary as a result of the following factors: 

1) rw, the nominal radius of the well (Q/s increases with 
increasing rw, i.e. directly proportional to log rw^) 

2) t, the pumping period (Q/s decreases with increasing 
time, i.e. inversely proportional to log t) 

3) sw, the well loss (Q/s decreases with an increase in 
punning rate in wells with a high sw) . 

The relationship between the theoretical specific 
capacity and the coefficient of transmissibility of an artesian 
aquifer is shown in Figure 5 (after Csallany et al, 1963) . It 
is assumed that the conditions depicted in Figure 5 are 
applicable in the study area. Permeabilities were computed 
from specific-capacity data, using Figure 5, for artesian wells 
in the study area and the results are presented in Table 9 
below. According to Walton (1962) and Walton et al (1962) , the 
coefficient of transmissibility, computed from specific-capacity 
data, does not vary significantly between a pumping period of one 
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TABLE 9. COMPARISON OF PERMEABILITIES ("K") DERIVED FROM SPECIFIC- CAPACITY AND PUMPING- 

TEST DATA OF THE SAME V'fELLS IN THE IFYGL DUFFINS CREEK STUDY AREA 



Material 



No. of wells Mean diameter 



(inches) 



Mean puinping 
period 

(hours) 



Mean "K" from 
specific 
capacity data 
(Igpd/ft2) 



Mean "K" from the 
modified non-equilibrium 
well formula- 
(Igpd/ft^) 



gravel 31 

sand & gravel 10 

sand 8 

till, clay, silt 4 

bedrock 42 



5.81 (1.18) 
5.65 (0.82) 
5.21 (0.68) 
5.87 (0.59) 

5.82 (1.04) 



4.06 (2.34) 

3.50 (1.58) 

3,12 (1.55) 

2.81 (1.47) 

3.23 (2.04) 



1357 (2118) 

842 (920) 

244 (210) 

240 (10 8) 

28 (39) 



4 38 (772) 

304 (44 3) 

73 (6 8) 

25 (29) 

7 (9) 



(1.18) standard deviation 



to 8 hours and a limited range of the nominal radius (rw) 
of the well. In order to minimize these errors, only wells 
in the study area with a pumping period ranging from 1 to 10 
hours and with a well diameter of 4 to 8 inches were used in 
the computations of permeability presented in Table 9. 

The data in Table 9 indicate that the permeabilities 
derived from the specific-capacity data are larger by 
approximately a factor of three than those permeabilities 
derived from the modified non-equilibrium well formula (Table 8) . 
Sufficient data on well losses (sw) and hydraulic properties 
of the confining beds are not available to delineate the 
permeability from the specific-capacity data to any greater 
degree of accuracy. 

GROUND-WATER CONTRIBUTION TO LAKE ONTARIO 

Introduction 

The Darcy formula, Q = T I L , was used to compute the 
amount of direct ground-water inflow to Lake Ontario, where: 

Q = quantity of water in Imperial gallons 

per day (IGPD) 
T = coefficient of transmissibility (I opd/ ft) 
I = hydraulic gradient (feet per foot) 
L = width of flow cross section (feet) . 

In order to evaluate the values to be used for the parameters 
listed above, several assumptions on the ground-water flow 
field were required. The assumptions used in the determination 
of ground-water inflow to Lake Ontario are listed below with 
respect to contributing area, depth of active flow, saturated 
thickness, transmissibility and hydraulic gradient. 

Contributing Area 

As indicated in an earlier section, the area near the 
Lake Ontario shore is considered to be in a zone of general 
ground-water discharge. In order to facilitate the computation 
of ground-water inflow to the Lake, it was assumed that the area 
within a distance of one mile from the shore was representative 
of the conditions in the zone of ground-water discharge. 

Depth of Active Flow 

The mean value of bedrock penetration for 112 wells 
completed in the bedrock within the study area was found to be 30 feet. 
Wells completed within the assumed zone of ground-water discharcre 
(contributing area) indicate a mean value of 10 feet for bedrock 



penetration. From these data it was assumed that the mean value 
would reflect the practical limit of ground-water exploration 
for domestic supplies and therefore, significant ground-water flow 
in the bedrock would not extend below 50 feet of bedrock pene- 
tration. This arbitrary value of 50 feet of bedrock penetration 
IS considered to be the boundary between the active and stagnant 
parts of the ground-water reservoir. 

Saturated Thickness 

The static water level reported in the water-well records was 
used to determine the saturated thicknesses of materials in wells 
within one mile distance of the Lake Ontario shore (Table 10) . In 
some wells, a relatively low, static water level in conjunction 
with a thin overburden thickness indicated that the contri- 
bution of ground water from the overburden was negligible. 

Transmissibility 

Permeability values , derived from short-term pumping 
tests (Table 8) using the modified non-equilibrium well 
formula (Cooper & Jacob, 1964) , were used to compute trans - 
missibilities in the zone of ground-water discharge. The 
various materials reported in the logs of wells within a 
distance of one mile from the shore were assigned a permeability 
which was then multiplied by the saturated thickness to deter- 
mine the transmissibility. The total transmissibilities for 
each well were summed and an average transmissibility of 930 
Igpd/ft. (Table 10) was determined for the materials within 
the zone of groiand-water discharge. 

Hydraulic Gradient 

The average hydraulic gradient, I, was assumed to be 
0.010 (Map 10) for the study area as indicated in an earlier 
section on ground-water movement. 

Computation of Ground-Water Inflow 

The values for T and I, as derived above, were siobstituted 
in the Darcy formula to obtain a value of ground-water con- 
tribution per mile length of shore line as follows: 

Q = TIL 

= 9 33 X 0.010 X 52 80 

= 49,262 Igpd per mile length 

= 0.092 cubic feet per second (cfs) per mile length. 

It is assumed that this value of ground-water discharge for the 
study area is representative of 35 miles of Lake Ontario shore- 
line from the mouth of the Humber River, near the western 
limits of the City of Toronto, to the City of Oshawa, giving 
a grand total of 3.2 cfs for direct ground-water discharge to 
Lake Ontario. Itiis value is similar to the values of 2.61 and 
3.29 cfs for roughly the same area as computed by Haefeli (19 72) . 



TABLE 10. MEAN VALUES OF TRANSMISSIBILITY (T) ,, STATIC WATER- 
LEVEL DEPTH, OVERBURDEN THICKNESS AND BEDROCK 
PENETRATION FOR WELLS COMPLETED WITHIN ONE MILE OF 
THE LAKE ONTARIO SHORE IN THE IFYGL DUFFINS CREEK 
STUDY AREA. 



Mean Over- Mean Bed- Mean 

burden roek Pen- Static 

NO'. of Mean T thickness etration Level 
Wells (gpd/ft.) (ft.) (ft.1 (ft.) 

wells completed 

in the 

overburden ti 1036(810) 33(14) ^ IS (11) 

wells completed 

in the bedrock 18 610(219) 45(22) 10*(14) 27*(21) 

all wells 74 933(735) 36(17) 18*(15) 

(219) standard deviation 

* anomalous well with 19 3 feet of bedrock penetration and 

static water-level of 234 feet was not considered. 

Additional Ground-Water Inflow 

Relatively high bluffs occur along the Lake Ontario 
shore, west of the study area, in the Borough of Scarborough. 
The height of tliese bluffs averages approximately 410 feet asl 
or 160 feet above lake level, along a nine-mile stretch of 
shoreline. Assuming that the discharge of ground-water may be 
occurring in the lower one-third (Meyboom et al , 1966) of these 
bluffs (53 feet) thie value of 3.2 cfs would then be minimal 
for the ground-water discharge to Lake Ontario from this hydro- 
geological region. 

A calculation of the additional discharge from this 
stretch of the Scarborough Bluffs was made using values of 
transmissibility derived from the Duffins Creek study. The 
average permeability per foot of saturated overburden thick- 
ness is 57.5 gpd/ft. (Table 10, "wells completed in the over- 
burden" column: 1036 -^ (33-15) = 57.5). The mean saturated 
thickness of the overburden in the study area is 18 feet 
(Table 10, "all wells" column: 36-18 = 18). The additional 
saturated overburden thickness in the Scarborough Bluffs 
section would then be 35 feet (subtracting the mean saturated 
thickness of the study area from the estimated value of 
53 feet of saturated thickness in the Scarborough Bluffs). 
Substituting these values in the formula: 

Q = TIL = KmIL = (57.5) x 35 x (0.010) x (5280 x 9) 

^ 954,700 IGPD = 1.77 cfs. 

Therefore, an additional 1.8 cfs of ground water may be dis- 
charging directly to Lake Ontario from the highest part of 
the Scarborough Bluffs giving a grand total of 5.0 cfs of 
ground-water discharge to Lake Ontario. 



SUMMARY AND CONCLUSIONS 

A hydrogeQlogical evaluation of the lower part of the 
Duff ins Creek drainage basin was made for the IFYGL program, 
in order to estimate the amomit of ground-water inflow to 
Lake Ontario, in the hydrogeologic region along the north shore 
of Lake Ontario from the Humber River near the western limits 
of the City of Toronto to the City of Oshawa (Figure 1) . The 
hydrogeology of the study area was assessed utilizing field 
investigations of the geology, test drilling to provide infor- 
mation on aquifer charaGteristics , subsurface geology and 
ground-water levels and information from water-well records on 
file with the Ministry of the Environment. 

Ifie study area lies in parts of the physiographic regions 
of the Iroquois Plain and South Slope (Map 1) . The Iroquois 
Plain is the lowland which borders Lake Ontario, that was 
inundated and smoothed by wave action of glacial Lake Iroquois . 
The South Slope region borders the study area to^ the northwest 
and consists of drumlinized till plain. Regional topographic 
gradients range from 25 feet per mile in the Iroquois Plain to 
150 feet per mile in the South Slope region. 

The climate of the study area is influenced by the 
presence of Lake Ontario to the south, A long-term annual 
precipitation and temperature of approximately 32 inches and 
44OF, respectively, is considered to be representative over 
the study area. 

The bedrock of the area is comprised of Upper Ordovician 
shales of the Whitby Formation which are overlain by Wisconsinan 
drift of the Pleistocene Epoch. Two major bedrock valleys 
were delineated in the study area and they are roughly coin- 
cident with present-day drainage. Subsurface information from 
water-well records on file with the Ministry of th.e Environment 
indicates that the area was inundated by high-level glacial 
lakes during Wisconsinan time. The permeable materials (sand 
and gravel sediments) of the Thorncliffe Formation that were 
deposited in these lakes are presently being utilized as a 
source of water for domestic and agricultural uses. Correlation 
of the unconsolidated materials in the study area with 
materials to the east (Bowmanville area) and to the west 
(Scarborough area) is given in Table 2. 

An excess of 4.29 inches of precipitation over the 
long-term annual average during the IFYGL period (April 1, 19 72 
to March 31, 19 73) resulted in a rise of the static level of 
0.85 feet in the ground-water reservoir of the study area 
(Table 6) . 

Water-level data from water wells drilled in the study 
area suggest that a hydraulic connection exists between the 
overburden materials and the upper part of the bedrock. The 
mean static water-level for overburden wells in the study area 
is 15.5 feet below the ground surface, and it is suggested 
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that the potentiometric surface conforms to the topographic 
expression of the land surface (Map 10) . Ground-water divides 
are generally coincident with surface-water divides indicating 
that ground-water movement is from topographically high to 
topographically low areas. Piezometric data indicate that the 
area near the Lake Ontario shore appears to be in a zone of 
general ground-water discharge whereas areas inland from the 
Lake appear to be in a general area of recharge (Table 7) , 

Secondary porosities developed in the upper part of the 
bedrock likely constitute the major source for ground-water 
storage in the bedrock. The contributing thickness of bedrock 
to significant ground-water flow was estimated to be 50 feet . 
The extent of granular deposits in the overburden, which are the 
most productive materials for ground water, were delineated 
for the study area (Maps 6 to 9 , inclusive). 

Permeabilities of the overburden materials and the bedrock 
were computed using the modified, non-equilibrium well formula 
and specific-capacity data (Table 9). Permeabilities derived 
from the specific- capacity data are larger by approximately a 
factor of three than those permeabilities computed from the 
modified, non-equilibrium well formula. Water wells in the 
overburden are more productive than wells completed in the bedrock 
(Figure 4) with overburden wells having a specific-capacity range 
of 0.8 to 5.0 Igpm per foot of drawdown and bedrock wells a 
specific-capacity range of 0.07 to 0.7 Igpm per foot of drawdown. 

Using permeability values derived from short-term pumping 
tests, an average transmissibility of 930 Igpd per foot was 
computed for materials within one mile of tlie present Lake Ontario 
shore in the study area. The mean overburden thickness penetrated 
by these wells was found to be 36 feet. Ground-water inflow 
to Lake Ontario, within the study area, was estimated to be 
0.092 cfs per mile length of shoreline. A basic value of 3.2 
cfs was estimated for direct ground-water discharge to 
Lake Ontario for the hydrogeologic region extending from the 
Hurriber River to the City of Oshawa along the north shore of 
Lake Ontario. An additional 1.8 cfs of ground-water discharge 
to Lake Ontario may be occurring along a 9-mile stretch of the 
comparatively high Scarborough Bluffs, west of the study area 
in the Borough of Scarborough. 
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